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ABSTRACT 

An investigation on correlations between photospheric current helicity and subsur- 
face kinetic helicity is carried out by analyzing vector magnetograms and subsurface 
velocities for two rapidly developing active regions. The vector magnetograms are 
from the SDO/HMl (Solar Dynamics Observatory I Helioseismic and Magnetic Im- 
ager) observed Stokes parameters, and the subsurface velocity is from time-distance 
data-analysis pipeline using HMI Dopplergrams. Over a span of several days, the evo- 
lution of the weighted current helicity shows a tendency similar to that of the weighted 
subsurface kinetic helicity, attaining a correlation coefficient above 0.60 for both ac- 
tive regions. Additionally, there seems to be a phase lag between the evolutions of the 
unweighted current and subsurface kinetic helicities for one of the active regions. The 
good correlation between these two helicities indicate that there is some intrinsic con- 
nection between the interior dynamics and photospheric magnetic twistedness inside 
active regions, which may help to interpret the well-known hemispheric preponder- 
ance of current-helicity distribution. 

Subject headings: Sun: photosphere — sunspots — Sun: helioseismology 



1. Introduction 



Parker] (I1955|) initiated a model for the cycle of magn etic field in the sun with a mecha- 
nism named "cyclonic motion", and ISteenbeck et al. Jl966h established the feasibility of AC or 



DC dynamos in astrophysical environment. At that time, they formulated a c onsistent expla 
nation of the role of the cy clonic motions, which was later called a-effect (e.g.. [Roberts & Stix 
197 IPouquet et al.1 (|1976|) suggested that turbulence and magnetic fields were common features 
of many celestial bodies and the a-effect contained both the kinetic and magnetic helicity. The 
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corre l ation b etween these tw o kinds of helicities has been under debate (e.g.. iBrandenburg et al 
1990. Il995|) . iKeinigsl (I1983|) gave an example of a-effect including zero kinetic helicity and the 
a-effect relied on the ensemble average of fluctuating velocity and magnetic field. That is to 
say, the generation of magnetic field requires current helicity but not necessarily the kinetic he- 
licity. The a-effect is central to dynamo theory and magneti c field generati on, however, it lacked 
observational and experimental supports over a long term (ISokoloffl 120071). Not until the recent 



20 years, mirror asymmetries of magnetic field (magnetic helicity) flS eehafer 1990; Pevtsov et al. 



19951 : iBao&Zhand 1 1 998: 



Hagino & Sakurai 2004 



Zhang et al.ll2010|) and the velocity (kinetic 



helicity) (|Zhao & Kosovichevl 120031) were found in observations. Apparently, frozen flow with 
the magnetic field in the fluid of high magnetic Reynold number is a necessary condition for the 
ct-effect in the framework of Parker's dynamo model. It can be easily inferred that the variation 
of kinetic and current helicity in the fluid might be a good manifestation of the freezing process. 
Although the speculation of the depth of a-effect is still under debate, it is expected that some 
relevant characteristic might be captured even near the photosphere. 

Local helioseismology provides a uniq u e tool to determine the sub-photospheric flows of 
active regions. A statistical study by Izhaol (|2004l) showed that the subsurface kinetic helicity 
inside active regions observed by the Solar Heliospheric Observatory I Michelson Doppler Im- 
ager (SOHO/MDl) observations seemed to ha ve a hemispheric preponderance, like wh a t magnetic 
(or cu rrent) helicity observations had shown ( Pevtsov et al. 1995 : Bao & Zhanglll998|) . Gao et al 



(|2009r) analyzed the connection between the photospheric current helicity, calculated from vector 
magnetograms observed by Huairou Solar Observing Station, and the subsurface kinetic helicity 
measured from MDI observations in 38 solar active regions. Although there was an opposite hemi- 
spheric trend between the sign of current helicity and that of subsurface kinetic helicity near the 
solar surface, the result did not support that the subsurface kinetic helicity had a cause and effect 
relat ion with the phot ospheric current helicity at the depth of 0-12 Mm. Similar result was reported 
by Ifvlaurya et all (1201 1|) as well. 



Now, Solar Dynam ics Observatory I Helioseismic and Magnetic Imager (SDO/HML lScherrer et al, 
2012l : ISchou etall2012h observations provide an unprecedented opportunity to investigate the con- 
nection between subsurface kinetic helicity and current helicity, as both subsurface flow velocity 
and photospheric vector magnetic field are available at the same time. Subsurface kinetic helic- 
ity can be computed from subsurfac e flow velocities, which are routinely processed through the 
HMI time-distance analysis pipeline (|zhao et al. 2012 ). Photospheric current h elicity is able to be 
computed from the photospheric vector magnetograms ( Hoeksema et al. 2012h . In this Letter, we 
focus our study on active regions NOAA AR 1 1 158 and AR 1 1283, both of which have continuous 
observational coverage that allows us to investigate the relationship between subsurface kinetic 
helicity and current h elicity over a lo nger period rather than just a snapshot comparison as done 
in our previous study (|Gao et alj|2009l) . In §2, we describe the procedure of observation and data 
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preparation for this analysis, and present our results in §3. We then summarize and discuss our 
results in 84. 



2. Observation and Data Reduction 



HMI observes the full-disk Sun continuously, providing Doppler velocity and line-of-sight 
magnetic field, a mong others, with a 45-sec cadence, and also vector magnetic field with a ca- 
dence of 12 min (|Schou et al.ll2012|) . Each full-disk image has 4096 x 4096 pixels with a spatial 
resolution of 0.504 arcsec pixel -1 (i.e., approximately, 0.03 heliographic degree pixel -1 at the so- 
lar disk center). Subsurface flow velocities are computed f rom HMI Doppl er-shift observations 
using time-distance helioseismology data-analysis pipeline (IZhao et al.ll2012l) . In practice, to exe- 
cute this data-analysis pipeline, users provide the Carrington coordinate for the interested area and 
time of the interested period, then the pipeline code selects an area of roughly 30° x 30° centered 
at the given coordinate, and a duration of 8 hours with the given time as the middle point. Nor- 
mally, the pipeline generates a subsurface flow field consisting 256 x 256 pixels with a horizontal 
spatial sampling of 0.12° pixel -1 , and a number of depths covering from the photosphere to 20 
Mm in depth. We performed analysis in all depths, and found that the shallowest depth, i.e., - 
1 Mm, gave the most consistent and robust results. Thus, in this Letter, we only present results 
obtained from this depth and leave analyses of deeper layers for future studies. And also, a recent 
comparison of subsurface flow field at the depth of - 1 Mm, obtained from this time- distance 
data-a nalysis pipeline, and photospheric flows obtained from the DAVE4VM techni que (Schunkl 
2008) found a reasonable agreement in horizontal flow fields inside active regions ( Liu & Zhaol 
20121) . Vector magnetograms are also from HMI o bservations, and these data have also been 
widely used in carrying out quantitative analysis (e.g. JSun et al.ll2012|) . 



NOAA AR 1 1 158 was a rapidly developing active region, in which an X 2.2 solar flare started 
at 01:44 UT and peaked at 01:56 UT of 2011 February 15. This active region, a s well as the X- 
class solar flare that occurred in it, were widely studied for different purposes (e.g., Sun et al.ll2012 : 
Jing et alJl2012clLiu & Zhadl2012|) . Figure [Ik shows an example of velocity field of AR 1 1 158 at 
the depth of - 1 Mm, with the background showing the vertical component of the flow and 
arrows showing the horizontal components. Figure [Tp shows an example of vector magnetogram 
of NOAA 11158 taken at a time similar to the subsurface flow field. Figure [2] shows the evolution 
of the averaged absolute values of three components of the magnetic field and three components of 
the subsurface velocity, respectively. It can be found that the magnitudes of the components of the 
velocity and magnetic field remain relatively stable, and there is no coherent or sudden change for 
any of these quantities. Meanwhile, the fluctuations of these quantities are also small, as shown in 
the curves of standard deviations. 
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Another active region, NOAA AR 1 1283, is also studied. This region was very flare produc- 
tive, generating two C-class solar flares from 22:00 UT to 23:33 UT on 201 1 September 5, an X2.1 
flare that started at 22: 12 UT and peaked at 22:20 UT on September 6, an XI .8 flare that started at 
22:32 UT and peaked at 22:38 UT, and an M6.7 flare on September 8 when the region was close 
to the west limb. As shown in §3, the same analysis on evolutions of photospheric current helicity 
and subsurface kinetic helicity is applied on both AR 1 1 158 and AR 1 1283, and similar results are 
achieved for both active regions. 



3. Results 

The main purpose of this work is to study whether there is a correlation between the pho- 
tospheric current helicity and subsurface kinetic helicity insid e active regions. Curr ent helicity is 



defined as: H c = B-(VxB). Similar to what was employed by lBao & Zhangl(ll998|) . the averaged 
value of the vertical component density of the current helicity, denoted as (if*), was used in this 
study. In addition, we also compute the weighted (if*), defined as (H*) divided by the total mag- 
netic field strength, i.e., (if*) /|B 2 |. This parameter has the same dimension as the mean twist a, 
a commonly used proxy of magnetic helicity. However, a is defined as the curl of transverse mag- 
netic field dividing the corresponding longitudinal field, hence carries a same sign as the weighted 
(if*)- Similarly, kinetic helicity is defined as if/. = v • (V x v), and we take the averaged value of 
its vertical component, (iff), and the value weighted by the speed, (iff)/|v 2 |, for further studies. 
To compute both values of (if*) and (iff), we only use the areas whe r e \B Z \ > 50 Gs, and this is 



a popular approach taken in many previous studies, e.g.. lBao & Zhan g (1998). The magnetograms 
are rebinned to the same resolution of the velocity map, and all quantities are averaged over the 
same spatial areas. In this Letter, we study the evolutionary relationship between (Hf) and (iff) 
inside the two flare-productive active regions introduced above. 

The top panel of Figure [3] shows evolutions of the weighted current-helicity density and 
weighted kinetic-helicity density obtained at a depth of - 1 Mm for AR 1 1 158. The time series 
from 2011 February 13 to 17 are analyzed. The weighted current-helicity density with a 12-min 
cadence is averaged over a 4-hr period in order to compare with the weighted kinetic-helicity den- 
sity, which is computed from the subsurface velocity obtained from an 8-hr data sequence with a 
4-hr time step. Error bars are plotted for the current helicity with corresponding standard devia- 
tions. The two helicity curves show very similar varying tendencies and the correlation coefficient 
attains 0.67. Furthermore, we separate the data series into two sections. For the decreasing phase 
before 06:00UT of 201 1 February 14, the correlation coefficient between the kinetic-helicity den- 
sity and the current-helicity density is as high as 0.84. But for the rising phase after 06:00UT of 
2011 February 14, the correlation coefficient drops to 0.52. The decrease of the correlation may 
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be related to the X2.2 flare event that occurred at 01:44 UT of 201 1 February 15. For the calcula- 
tion of these correlation coefficients, the included data points n (correlation coefficients r) are 22 
(0.67), 7 (0.84), and 15 (0.52), thus the degrees of freedom are 20, 5, and 13, respectively. The cor- 
responding critical values of the Pearson correlation coefficients for a significance level of 0.95 are 
0.423, 0.754, and 0.514, respectively. Therefore, the correlation coefficients from our calculations 
are all highly significant. 

The bottom panel of Figure [3] shows the evolutions of the unweighted current-helicity density 
and unweighted subsurface kinetic-helicity density. Different from the similar evolutionary ten- 
dency of the weighted parameters, the two unweighted parameters seem to evolve out of phase. As 
shown in this panel, the unweighted kinetic helicity seems to have a decreasing phase about 8 hr 
earlier than that of the current helicity before 10:00 UT of February 14, and have an increasing 
phase about 4 hr behind that of the current helicity after 18:00 UT of February 14. The bottom 
panel of Figure [2] shows the evolution of B 2 and v 2 . These two quantities seem to evolve out of 
phase at the beginning and more in phase after mid-day of February 14. This may account for the 
evolutionary difference in the weighted and unweighted kinetic helicty. 

The snapshots in Figure |4] show maps of the weighted and unweighted (H%) as well as the 
weighted and unweighted (H') for some selected periods. It can be found that both of the kinetic- 
helicity parameters are more fragmented than the current-helicity parameters, but there is not clear 
correspondence between the sign distributions of the kinetic helicity and the current helicity. Be- 
sides, we also plot weighted (Hf) and unweighted {H*) after binning down to match the spatial 
resolution of (H£). Furthermore, we compute the correlation between the maps of weighted (un- 
weighted) current and kinetic helicity with the same resolution. The values are 0.019 and 0.028 
respectively. 

For the other active region NOAA AR 11283, the evolutions of both current helicity and 
subsurface kinetic helicity, both weighted and unweighted, are shown in Figure |5] Similar to the 
results for AR1 1 158, the evolutionary curves of the weighted parameters have a better correlation, 
which is 0.62 for this case, than those of the unweighted parameters. However, for the unweighted 
parameters, we do not see the phenomenon of that the subsurface kinetic helicity evolves ahead 
of the current helicity before the flare and behind the current helicity after the flare, like what is 
observed in AR 1 1 158. The evolution of B 2 and v 2 for this active region does not show any phase 
lag either. 
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4. Discussion 

We have studied the correlation between the evolutionary curves of the current-helicity den- 
sity, computed from HMI vector magnetic field, and the subsurface kinetic-helicity density, com- 
puted from subsurface velocity field obtained from HMI time-distance data-analysis pipeline, for 
two flare-productive active regions. For the weighted helicities, we find that the evolution of the 
current helicity and subsurface kinetic helicity has a high correlation, larger than 0.60 for both ac- 
tive regions. But for the unweighted case, essentially, both helicities do not show a high correlation 
for the studied period for both regions. It is not clear why the weighted parameters show a high 
correlation while the unweighted parameters do not. Still, it is very interesting to see the weighted 
helicities attain a high correlation during the evolution of both active regions, especially when 
considering that these two helicities are calculated from very different data. Although the original 
observations were both from the same instrument HMI, current helicity is calculated from vector 
magnetic field which was derived from the 12-min-cadence Stokes parameters, and the subsurface 
kinetic helicity is calculated from the 45-sec-cadence Dopplergrams after some very complicated 
helioseismological processing. The high correlation between these two helicities indicate there 
are indeed some intrinsic connections between the subsurface kinetic helicity and the photospheric 
current helicity. 

Since the hemispheric prepondera nce of magnetic (or current) helicity was reported (e.g., 
Pevtsov et al.lll995 : Bao & Zhanglll998 ). many efforts have been made to explain this observational 



phenomeno n based on the solar dynamo (e.g.. lKleeorinl l2003 



Choudhuri et al 



2006 . 2012 ). and the S-effect of emerging magnetic flux (|Longcope et a~ 1998 ). The adopted 



2004: IZhang etal. 



ve- 
locity field here mainly comes from the region near the solar surface, so the correlation between 
the short time scale variations of kinetic and current helicity may be suitable to be explained by 
using the S-effect. The S-effect suggested that during the rising of the magnetic flux from convec- 
tion zone to the photosphere, the large-scale writhing and the local-scale twisting together shape 
the tilt angle and magnetic (or current) helicity distributions of active regions. The subsurface ki- 
netic helicity used in our analysis is presumably corresponding to the twisting inside active regions 
in the upper convection zone. The in-phase evolution of the subsurface kinetic helicity and the 
current helicity probably indicate that the twisting beneath active region surface indeed play an 
important role to shape the current helicity distribution observed in the photosphere. On the other 
hand, perhaps it is not surprising that these two helicities evolve in phase, because beneath the pho- 
tosphere magnetic field lines are frozen with plasma and the evolution of photospheric magnetic 
field somehow reflect the subsurface motions. 

Another fact we cannot ignore is that despite the high positive correlation in the evolutionary 
curves of the two helicities, the signs of the two helicities do not often stay the same. In particu- 
lar, for AR 11158, the signs of the two helicities are more often opposite than same. Therefore, 
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although the evolution of the two helicities is in positive correlation, it is quite likely that we get 
a negative or n o correlation if w e pick just one random snapshot. This may help explain why our 
previous study (|Gao et al.ll2009() . which statistically analyzed the correlation of the two helicities 
using snapshots from many different active regions, did not find a correlation. The results presented 
in this paper also give us some guidance of how a statistical study with many active regions should 
be carried out, and we plan to perform such a statistical study once HMI accumulates sufficient 
number of active regions. 



Moll et al.l (|2012l) compared the flow structures in the upper photospheric layers from their 
simulations and found that vertically extended vortices were associated with magnetic flux con- 
centrations. Small-scale vortices occurring in simul ations of non-magneti c or weakly magnetized 
regions of the lower atmosphere was also studied (|Kitiashivili et al.ll2012|) . However, these simu- 
lations focused on regions outside of sunspots with weaker magnetic field strength, and the spatial 
scale was also smaller than the observational resolution. Thus, we believe our study provides a 
useful window to understand the temporal and spatial interactions between magnetic field and flow 
field in a larger scale, which can be hardly achieved using numerical simulations for the time being. 

It is certainly of great importance to study whether there is any indication of solar flare oc- 
currences hi&eeyohitwn^y curves of both the current and subsurface kinetic helicities (e.g., 
Bao et al.ll 19991 ; iKomm et al.ll2005|) . From the limited number of active regions analyzed in this 
study, it is premature to draw a conclusion. Figure[3b seems to show that for AR 1 1 158, before the 
flare occurrence the subsurface kinetic helicity evolved ahead of the current helicity, and evolved 
after the current helicity after the flare. However, this phenomenon is not observed in AR 1 1283. 
Nevertheless, this also needs a statistical study with more active regions and more solar flares. 
The availability of simultaneous HMI vector magnetic field and subsurface flow fields for active 
regions makes such a study possible. 
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Fig. 1. — (a) Example of subsurface velocity field of NOAA AR 11158 at a depth of - 1 Mm. 
Background shows the vertical component of the velocity and arrows show the horizontal compo- 
nents. The maximum vertical velocity is 285 m s _1 and the longest arrow represents a horizontal 
speed of 441.8 m s _1 . The spatial sampling of the vertical component is 2.016"x2.016", and the 
field of view is 3.02'x2.02'. (b) Example of vector magnetogram of NOAA AR 1 1 158, with white 
showing positive polarity and black showing negative. The longest arrow represents a transverse 
field of 2356 G, and the transverse field lower than 300 G is not shown. 
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Evolution of Magnetic Field 
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Fig. 2. — Temporal evolution of the three components of the magnetic field (top), the three com- 
ponents of the subsurface velocity field (middle) and the square of magnetic field and velocity 
(bottom). 



-12- 




Fig. 3. — Top: Evolution of the weighted current helicity (marked with star, corresponding to the 
left vertical axis) and the weighted subsurface kinetic helicity (marked with cross, corresponding 
to the right vertical axis) for AR 11158. Bottom: Same as the top panel but for the unweighted 
current helicity and unweighted subsurface kinetic helicity. Dashed lines, marked as "A" and "B", 
correspond to 10:00 UT and 14:00 UT of 201 1 February 14, respectively. 
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Fig. 4. — Snapshots showing weighted (iff) {upper left) and unweighted (iff) {upper right) at 
18:00 UT of 2011.02.15, weighted (ff 2 ) (Mz<Me and unweighted (ff 2 ) (AfwWZe rigfe) at 
19:00 UT of 2011.02.15 with original spatial resolution, and weighted (if 2 ) {Bottom left) and 
unweighted (if 2 ) {Bottom right) after binning down to match the spatial resolution of (iff). The 
parameters are displayed relative to the corresponding standard deviations, and color bars at the 
left lower corner of each panel show this scale. 
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Fig. 5.— Same as Figure H but for NOAA AR 1 1283. 



